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Abstract. We investigate the effects of disorder on driven lattice gases with open 
boundaries using the totally asymmetric simple exclusion process as a paradigmatic 
example. Disorder is realized by randomly distributed defect sites with reduced 
hopping rate. In contrast to equilibrium, even macroscopic quantities in disordered 
non-equilibrium systems depend sensitively on the defect sample. We study the current 
as function of the entry and exit rates and the realization of disorder and find that it is, 
in leading order, determined by the longest stretch of consecutive defect sites (single- 
bottleneck approximation, SBA). Using results from extreme value statistics the SEA 
allows to study ensembles with fixed defect density which gives accurate results, e.g. 
for the expectation value of the current. Corrections to SBA come from effective 
interactions of bottlenecks close to the longest one. Defects close to the boundaries 
can be described by effective boundary rates and lead to shifts of the phase transitions. 
Finally it is shown that the SBA also works for more complex models. As an example 
we discuss a model with internal states that has been proposed to describe transport 
of the kincsin KIFlA. 
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1. Introduction 

Driven diffusive systems play an important role in non-equilibrium statistical physics 
[H, [2] . Due to broken detailed balance they allow to investigate non-equilibrium effects. 
In addition, they serve as models for various transport processes ranging from vehicular 
traffic [31 H] to biological transport by motor proteins [3, EJ [7] . The paradigmatic model 
is the totally asymmetric simple exclusion process (TASEP) which was first introduced 
to describe protein polymerization in ribosomes [8]. It was solved exactly P, [10] which 
allows to study generic properties like the phase diagram and boundary induced phase 
transitions [11] that also occur in more complex driven systems without resorting to 
approximations or simulations. 

In contrast to the homogeneous case, much less is known for the TASEP with 
inhomogeneous hopping rates, i.e. disorder. Here in principle one has to distinguish 
particle-dependent and site-dependent hopping rates. The former case is simpler since it 
can be mapped onto an exactly solvable zero-range process [121 lEl HH fisjll . Therefore, 
at least for periodic boundary conditions, this case is well understood. In contrast, 
in the latter case exact results are known only for a single defect in an otherwise 
deterministic system with sublattice-parallel dynamics [T71 IT8l IT9] . For the general 
case, several numerical and (approximate) analytical investigations [20l [211 1221 [23] 
have revealed interesting behaviour already for single defects and periodic boundary 
conditions. However, far less is known for finite defect densities or open boundary, and, 
especially a combination of both. 

Apart from the fundamental theoretical interest in disorder effects in nonequilibrium 
systems, which are far from being well understood [24J, these are also of direct relevance 
for applications. Typical examples are found in intracellular transport processes where 
molecular motors often move along heterogeneous tracks like DNA or mRNA (see e.g. 
[25]). 

In this work we consider systems with binary disorder, i.e. with transition rates 
that can take two possible values that are randomly assigned to the sites. Sites with a 
lower transition rate are called defect sites (or slow sites due to their influence on the 
average velocity), while those with the higher transition rate are non-defect sites. 

For periodic boundary conditions the influence of single defect sites has been 
clarified by Janowsky and Lebowitz [211 [22]. The current through the system is limited 
by the transport capacity of the defect site leading to a density-independent current 
at intermediate densities. The corresponding steady state phase-separates into a high 
and low-density region separated by a shock. In contrast, at high and low densities the 
system is only affected locally by the presence of the defect. 

The effects of finite defect density in a periodic system have been studied in detail 
in [26l [271 [23] , mainly numerically. For binary distributions of defects, Tripathy and 
Barma [26l [27] have classified two different regimes: 1) a homogeneous regime with a 
single macroscopic density and non-vanishing current, 2) a segregated-density regime, 

X For investigations of disordered zero-range processes, see e.g. [16) . 
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with two distinct values of density and non- vanishing current. Considering the partially 
ASEP where disorder was realized by inhomogeneous hopping bias, they also found a 
vanishing-current regime which shows two distinct densities, but with a current that 
vanishes asymptotically for {L oo). They argue that phase separation can be 
understood by a maximum current principle: For a given mean density the system 
settles in a state which maximizes the stationary current. Thus the largest stretch of 
slow bonds acts as current limiting segment. For the same system, Juhasz et al. [28] 
introduced an effective potential and determined trapping times in potential wells to 
investigate the vanishing of the current in a finite-size scaling. 

In the case of open boundary conditions, not only detailed balance but also 
translational invariance is broken. Although in principle one expects the same phases 
as in the pure system, the phase boundaries and the nature of the transitions might 
change. For equilibrium systems the Harris criterion [29] allows to decide whether 
critical behaviour is altered by weak disorder. For nonequilibrium systems no such 
general statements are currently available [21]. 

As in the periodic case, for a single defect site generically disorder-induced phase 
separation into macroscopic regions of different densities is observed in certain parameter 
regimes. The presence of defects leads to a decrease of the transport capacity and the 
maximum current phase is enlarged compared to the pure system [30]. These results 
have been generalized to systems with a single stretch of consecutive defects, called 
bottleneck, or two defect sites (or bottlenecks) [31], [32], [SS] [M]. It has been shown that 
for such systems reliable analytical approximations exist which go beyond a simple 
mean-field approach to take into account the relevant correlations [311 . 

For open systems with a single defect (or bottleneck) the current depends on the 
position of the defect [32l |33l [31] if it is located close to a boundary. This edge effect 
due to the interaction between the defect and the boundary also occurs in systems with 
many defects and can be accounted for by effective boundary rates [34j . 

The focus of most previous investigations was on individual realizations of defect 
distributions while statistical properties of defect ensembles were not considered. 
Numerical investigations on ensembles have been made in [35], where the influence 
of defects on the phase transition between low- and high-density phase of the randomly 
disordered TASEP with open boundary conditions was studied. It was shown that the 
position of this phase transition is sensitively sample-dependent, even for large systems. 
This effect is due to defects near the boundaries which is consistent with the results in 
[34] . Krug [23] conjectured that also the maximum current is sample-dependent and is 
mainly determined by the longest bottleneck. In [32l [34] this was shown to be correct, 
at least for two bottlenecks which are not too close to each other. These observations 
lead to the Single Bottleneck Approximation (SB A) which is supported in this work by 
numerical and analytical arguments. 

For applications to real systems, macroscopic parameters and quantities are most 
relevant. Since we have seen that macroscopic quantities can depend on microscopic 
ones that can differ for different defect samples, we are mainly interested in determining 
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statistical properties, e.g. probability distributions and expectation values, of relevant 
quantities taking an ensemble of systems rather than looking at single samples. We 
therefor consider in this work a large ensemble of individual finite but large systems, 
while the individual values of these quantities might vary for different samples. 

The goal of this work is to understand the phase diagrams of driven lattice gases and 
give quantitative approximations for the expectation values of critical parameter values 
and the maximum current. Using Monte Carlo (MC) simulations, we therefore check the 
validity of the SBA and the concept of effective boundary rates on individual samples 
not only in the disordered TASEP, but also in a more complex model with internal 



states (NOSC model without Langmuir kinetics [Ml ET], see Appendix A ), which is a 
model for intracellular transport with KIFIA motor protein^. With the help of extreme 
value statistics these principles are used to derive approximations for expectation values 
(Sec. Hj). After checking the accuracy of the SBA we discuss in Sec. |5]the relevance of 
various possible corrections, e.g. through edge effects or effective interactions between 
the bottlenecks. Finally, in Sec. [6] the influence of the disorder on the phase diagram is 
investigated in more detail. 

2. Model definition 

Although our considerations are rather general, we focus here on the prototypical driven 
lattice gas, the totally asymmetric simple exclusion process (TASEP). The TASEP is 
defined on a lattice of L sites which can either be empty or occupied by one particle. A 
particle at site j can move forward to its neighbour site j ' + 1 if this site is empty. The 
corresponding hopping rate is denoted by pj. At the boundary sites j = I and j = L 
particles can be inserted and removed, respectively. If site 1 is empty a particle will be 
inserted there with rate a. On the other hand, if site L is occupied this particle will 
be removed with rate (3. Here we will use a random-sequential update corresponding to 
continuous-time dynamics. 

In a schematic form the transition rules read 

Forward hopping in the bulk: 10 — * 01 with probability PiAt 
Entry at the left boundary: — * 1 with probability a At (1) 
Exit at the right boundary : 1 — > with probability j3At . 

As the indices indicate, the hopping rates Pi are site dependent in the disordered TASEP. 
Other transitions are prohibited. 

We will also investigate a generalization of the TASEP where each particle can 



be in two different states 1,2. This model, which is defined in Appendix A has been 
proposed to describe the dynamics of KIFlA motor proteins on microtubules [361 EZj- 
In the NOSC model the forward-rebinding rate Uf is a parameter controlling the average 
velocity of single particles for which we allow disorder. 

§ These motor proteins belong to the kinesin family. 
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In this paper we focus on binary disorder for which the hopping rates on each site 
can take the two possible values p and q < p that are randomly distributed by the rule 



The parameter is the defect density. Sites with reduced hopping rate q will be called 
defect sites, or slow sites, while those with hopping rate p are non-defect sites or fast 
sites. In the following we will set p = 1 which can always be achieved by rescaling time. 

For convencience a sequence of / consecutive defect sites will be called bottleneck of 
length I in the following. A bottleneck of length / = 1 corresponds thus to an isolated 
slow site. 

We will focus on "finite but large" systems here, i.e. we neglect terms of magnitude 
0{1/L), but keep terms (9(1/ In L). This is motivated by the facts that (a) the maximum 
current decreases with increasing bottleneck length and (b) the length of the longest 
bottleneck grows logarithmically in L. 

3. Single Bottleneck Approximation 

Besides the macroscopic structure of the stationary state in dependence of the system 
parameters, i.e. the phase diagram, the main focus of our investigation will be the 
(stationary) maximum current J* for fixed bulk parameters p and q: 



In analogy with the terminology used in traffic engineering, we will call J* the ( transport) 
capacity. Besides the macroscopic structure of the steady state, the transport capacity 
will serve as main quantitative indicator for disorder effects. Furthermore the critical 
values etc and jSc where the transport capacity is reached are of interest. In the pure 
TASEP one has J* = J{a = 1, (3 = 1) = 1/4 and = f3c = 1/2. Both will change in 
the presence of defects and, based on previous results, could even be sample dependent. 

Investigations in several works [231 ISSl [M] indicate that in the TASEP with many 
defects, the longest stretch of consecutive defects (bottleneck) is the quantity that 
contributes most to the transport capacity. This is plausible if one assumes a local 
character of the bottlenecks by characterizing them by an individual transport capacity 
J*(/) depending on the length / and (possibly) position j. In the stationary state 
the total current is constant in space and is restricted by all bottleneck capacities, 
i.e. it can not exceed the minimum of all JJ{1). Since the transport capacity is 
decreasing monotonically with bottleneck size as was shown in [M], the minimum of Jj*(/) 
corresponds to the transport capacity J* {I*) of the longest bottleneck which consists of /* 
consecutive defects. Smaller bottlenecks do not contribute much as long as they are not 
too close to the longest one [31] • This motivates the Single Bottleneck Approximation 




with prob. (j) 
with prob. 1 — (f) 



(2) 



J* = max J(a, /3) . 



(3) 



(SB A): 
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n n pp 


Ipri p"f h 




•^SBA 


"^ISA 


1000 


0.05 


2 


2 


1 


0.2174 


0.2294 


0.2229 


1000 


0.1 


3 


12 


1 


0.2080 


0.2131 


0.2080 


1000 


0.2 


3 


2 


2 


0.1963 


0.2131 


0.2080 


3000 


0.1 


3 


4 


1 


0.2048 


0.2131 


0.2084 


3000 


0.2 


5 


5 


1 


0.1866 


0.1925 


0.1901 



Table 1. Comparison of Monte Carlo (MC) and SB A results for the transport capacity 
J* in the disordered TASEP with different system sizes L and defect densities (j). The 
transport capacity was obtained by Monte Carlo simulations for a = P = 0.5 
(column 6) for fixed slow hopping rate q = 0.6. This is compared with MC results 
(Jgg^, column 7) and the results obtained by ISA [34](Jjg^, column 8) for a single- 
bottleneck system with one bottleneck in the bulk whose length is the same as the 
longest defect in the simulated disordered TASEP (column 3). Columns 4 and 5 give 
the distance and length of the bottleneck next to the longest one. 



L 





/* 


distance 


length 


T* 


•^SBA 


1000 


0.05 


2 


4 


1 


0.07923 


0.08179 


1000 


0.1 


3 


2 


1 


0.07451 


0.07643 


1000 


0.2 


6 


3 


1 


0.06659 


0.06717 


3000 


0.1 


4 


6 


1 


0.07205 


0.07213 


3000 


0.2 


6 


3 


1 


0.06677 


0.06717 



Table 2. Same as in Table [1] but for the NOSC model without Langmuir kinetics. 
The forward hopping rate is inhomogeneous with uj^^^At — 0.58 and uj^^°^At = 0.32. 
The other parameters are fixed: uihAt — 0.8, uOgAt = 0.22, ujt, = 0. 



The transport capacity J* of a disordered system with randomly distributed 
defects is the same as the transport capacity J^ba ^ system with a single 
bottleneck if the length of this bottleneck is the same as that of the longest 
bottleneck in the disordered system. 

A similar conjecture has been made by Krug for periodic systems |23j . 

The SBA reduces the problem to the much simpler one of a single bottleneck in a 
system. In particular for the TASEP, efficient methods have been developed recently, 
namely the finite segment mean field theory (FSMFT) [31j and the interaction subsystem 
approximation (ISA) [M] . 

We expect the SBA to work for generic driven lattice gases, especially for low defect 
density 0, where the average distance between defects is large and their interactions can 
be neglected. As an example we have tested it not only for the TASEP, but also the 
disordered NOSC model in the limit of vanishing Langmuir kinetics. In both systems 
the average velocity of the particles is dependent on one or more transition rates. In 
the TASEP the hopping rate p is such a parameter, while in the NOSC model the 
forward- rebinding rate cu/ is a parameter controlling the average velocity. 
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First we consider a fixed realization of disorder witli small defect density 0. In 
this case we have a system with dilutely distributed bottlenecks of different lengths. 
We want to test the SBA for the disordered TASEP and the NOSC model. For this 
purpose we simulated systems with different disorder samples and compared the results 
for the transport capacity J* with numerical and analytical results of systems with 
single bottlenecks in Table [H For each sample we identified the longest bottleneck 
/* and calculated the transport capactity ^sba(^*) ^ single-bottleneck system with 
just one bottleneck of size /*. One observes a quite good agreement, although the SBA 
seems to overestimate the transport capacity systematically. This is not surprising since 
effective interactions of the bottlenecks will lead to an additional decrease the current. 
From the results in |34j we expect that the main effect comes from bottlenecks near the 
longest one. There it was shown that for systems with two bottlenecks that, although 
the main reduction of the transport capacity comes from the longer one, the transport 
capacity will further be reduced if the distance between the bottlenecks is small. To 
illustrate this effect we included the distance of the nearest bottleneck in Table [H Since 
it is more probable to find a bottleneck close to the longest one for larger defect density 
0, the results tend to be less accurate with increasing 0. 

Surprisingly it seems that the values Jfg^ obtained by the semi-analytical ISA 
method [33] are more accurate than the numerical ones (J^ba) single-bottleneck 
system. This is because ISA usually underestimates the value of J*{1) in the TASEP 
with one bottleneck, while SBA overestimates the current. Thus errors cancel. 

4. Probability distributions and expectation values in SBA 

As we have seen, the transport capacity depends quite strongly on the particular sample 
of the defect distribution, i.e. the size of the longest bottleneck. Usually in real systems 
the exact distribution of defect sites is not known, particularly the size and position of 
the longest defect can not be identified. Then a statistical treatment, i.e. considering 
an ensemble of systems with fixed defect density, is more appropriate. It allows to 
determine expectation values for quantities like currents and effective boundary rates 
(see Sec. 16. ip . This is especially relevant for applications e.g. to intracellular transport 
on cell filaments. Each cell consists of a large number of filaments that serve as tracks 
for motor proteins, and often inhomogeneities play an important role [25]. Therefore 
each filament can be modeled by a driven lattice gas on a linear chain [361 EH] and the 
quantities of interest are averages rather than the properties of individual chains. 

In this section we want to approximate the expectation value of the transport 
capacity J*{q, 0, L) for fixed defect density and finite but large system size L. In the 
last section we have shown that for small the capacity depends approximately on the 
size of the longest defect. Therefore we first determine the expectation value for the size 
of the longest defect in such a system. 

We now consider a given sample at defect density and system size L. The k- 
th bottleneck has length Ik and in the following two consecutive fast sites j, j + I will 
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be interpreted as a bottleneck of length / = located at site j. This implies that the 
number A*";, of bottlenecks is equal to the number Nf of fast sites, since each bottleneck is 
followed by a fast site (JJ]. The bottleneck length / is a random variable with distribution 

P^(/)=0'(l-0). (4) 

Since on average the fraction of fast sites is (1 — 0), the mean number of fast sites is 
{Nf) = (1 — 4>)L. The length of the longest bottleneck is I* = ma.x{lk\k = l,...,Nf}. 
The statistics of the maximum of independently distributed random values is governed 
by extreme value statistics [39]. It says that for a continuous probability distribution 
P{1) that decays exponentionally or faster for I oo, the probability density of /* 
being the maximum value of independently distributed random values is for large N 
asymptotically described by the Gumhel distribution ^39j 

G{u) = e-^e-'^"" (5) 

where u = u{l*) is a rescaled and shifted function of /* depending on the details of the 
probability distribution P{1). 

However, since in our case the probabihty distribution is discrete we need to be 
careful. Therefore, following the derivation used in [39] for continuous distributions, we 
derive the probability distribution of the maximal bottleneck length explicitly in order 
to control errors made by approximations. This will also provide an explicit expression 
for u{l*). 

The probability of a bottleneck being shorter than /' is 

V-l 

Since the Ik are independently distributed, we have the probability that all Ik are smaller 
than /': 

H^{1') = P^ilYf = exp (iV^ ln(l - 0^) 

= expf(l-0)Lln(l-0O)) . (7) 



For large L this probability is significantly larger than zero only for -C 1 and we can 
use the approximation ln(l — 0'') ^ —0'', thus 

^exp(-0''(l-0)L). (8) 

As was shown in [39], the error of this correction is (9(1/L^) for exponential P{1)- Thus 
we can neglect finite size corrections. 

The probability that all values are smaller than I' is equal to the probability that 
the maximum /* is smaller than /', 

I'-i 

H,{i') = j2nn- (9) 



II We neglect the possible exception at the right boundary. 
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V{1*) is the probability that the longest bottleneck has length /* which is explicitly given 

by 

v{r) = H<{r + i)-H<{r) (lo) 
= H'^{r + ]^) + o {{M*f) 

^ - L(l - 0)0'*+^ ln0exp - 

= -ln0e-"e-""" = -ln0G(u) (11) 
where G{u) is the Gumbel distribution and we have introduced the function 



u{l*) = - + - J In - ln(l - 0) - In L . (12) 

Now we assume this probability distribution to be continuous. Using the Euler- 
Maclaurin formula the expectation value of /* becomes 

'dl* 



(D = Vrp(/*) ^ [ i*v{i*] 

I* -^0 



1, ^ 



j\^{r + l)di* (13) 



since H^{1* + 1/2) is the cumulative distribution function of V{1*). Therefore we have 
lim^.^oo -f^<(^*) = 1 and it is bounded. Hence 

(r) = Lij<(L + -)- / H^{i*)dr 

^ Jl/2 

= LH^L + 1/2) - ^-'-^'^ - ^'^'"*> - - (14) 

In 

where Ei(x) is the exponential integral function. It can be expanded ^0] : 

Ei(— x) = 7e + ln(a;) + 0{x) for small x > 
Ei(— x) = 0{e~^) for large |x|, x > 

where 7e ~ 0.5772 is the Euler-Mascheroni constant. 
With these expansions we have for large L 

7e + ln(L(l - 0)0^+5) + C(e-i'^<^^) - C(e-^(i-'^)v^*)) 



(15) 



(/*) = L- 



In. 



7e + In L + ln(l - 0) + In 0(L + 1/2) + 0{{) 

= ^ M 

For finite but large systems, we can neglect terms of order 0{1/ L). Thus we obtain for 
the expectation value of the longest bottleneck 

_ lnL + ln(l-0)+7e 1 

ln(l/0) 2- ^''^ 

{I*) diverges for infinite systems, as expected. However, it grows only of order O(lnL), 
so that we have to keep this term in finite but large systems. 
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L 





no. of samples 


{J*)mc 


(-^sba)mc 


500 


0.1 


200 


0.2099 


0.2244 


1000 


0.2 


100 


0.1918 


0.2024 


3000 


0.1 


100 


0.2018 


0.2110 


3000 


0.2 


50 


0.1866 


0.1960 



Table 3. Comparison of disorder averages in MC results and SB A results for the 
expectation value of the tranport capacity. The defect hopping rate is q = 0.6. 
Column 4 shows the numerical results from MC simulations of the disordered TASEP. 
Column 5 displays results by SB A using the probability distribution (fTT|) . 



L 





no. of samples 


{J*)mc 


('^sba)mc 


500 


0.1 


200 


0.07495 


0.08010 


1000 


0.2 


100 


0.06852 


0.07258 


3000 


0.1 


100 


0.07438 


0.075553 


3000 


0.2 


50 


0.06852 


0.07258 



Table 4. Same as in Table [3] but for the NOSC model without Langmuir kinetics. 
The forward hopping rate is inhomogeneous with w^'^'^Ai = 0.58 and uj^^°'^At = 0.32. 
The other parameters are fixed: uoh^t — 0.8, cj^At = 0.22, — 0. 



If we approximate the transport capacity J*(</') for small by the corresponding 
current JgBA ^ system with one bottleneck, the expectation value is given by 
(J*(0)) = X]i^=o "^sba(^*) '^(^*)- Due to the approximation by a continuous function, 
the norm X]/*'=o ■^(^*) ^^"^ significantly deviate from one. In order to reduce this 
error we divide the result by 'Ya*=q'^{^*) 

/ 7* _ I^;*=o'^sba(^*)'^(^*) /.on 

(-'sbaHw - — — v{F) ■ ^ ' 

We can now either take numerical values for J* (I*) or (semi-) analytical ones from [3l] 
or [31]. Since V{1*) decays fast around (/*) it is sufficient to take into account only few 
terms in (fTSll in the vicinity of (/*). 

In order to display the generic character of the SBA, we show results for the 
transport capacity not only for the TASEP but also for the disordered NOSC model 
without Langmuir kinetics (Tables [T] and [2]). We observe a good agreement in both 
systems while the errors are of the same magnitude as for individual samples. This 
indicates that the probability distribution function for the longest bottlenecks is an 
appropriate approximation. 

5. Corrections to SBA 



In the following we consider corrections to the SBA and check the quality of this 
approximation and the range of its validity by statistical means. 
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In principle, corrections to the transport capacity could come from the following 
effects: 

• The longest bottleneck (length /* = max{/i,/2, ■ ■ ■}) is located near the boundary, 
not in the bulk as assumed in SBA. Since the probability that a bottleneck at a given 
site is smaller than / is P<{1) = 1 — 0' (see ([6])), the probability of finding the first 
longest bottlenecli^ of length / at distance x from a boundary is P{x) = (1 — 0')^0'. 
Therefore the average distance of the longest bottleneck can be approximated as 

^ 0-<'*) = 01/2^(1 - (f))e^^ = 0{L) (19) 

where we approximated the longest bottleneck by its expectation value ffT7|) . 
That means for large systems the longest bottleneck is, on average, far from the 
boundaries. However we see that for finite systems and small defect densities 
-C 1, (x) is becoming small, so that the boundaries might affect the transport 
capacity. This is due to a kind of "degeneracy" of the longest bottleneck, since for 
small defect densities the probability that there are many longest bottlenecks is 
high (e.g. for I* = 1 this degeneracy is 0{L)), thus contributions of samples with 
a longest bottleneck near a boundary are relevant. We therefore expect deviations 
from the SBA for very small defect densities in finite systems. The effect should 
vanish in the limit L ^ oo for fixed 0. 

• Other smaller bottlenecks near the boundary can be treated by introducing effective 
boundary rates (see next section). 

• Corrections from other bulk defects, i.e. "defect-defect interactions". Candidates 
for the leading contribution from this type of correction would be a) other long 
defects, i.e. defects of length I < I*, and b) defects (of arbitrary) located in the 
neighbourhood of the longest one. The results of [3lj indicate that the second 
correction is more important. 

In order to estimate the corrections we consider an ensemble of systems which all 
have a longest bottleneck of length I* and defect density 0. The slow hopping rate q is 
considered to be fixed. The longest bottleneck (or one of them in case of degeneracy) 
is located at an arbitrary position and the distribution of the other defect sites is not 
restricted. For this ensemble the average transport capacity is given by 

(j*)(0,L,r) = 5^'j;(x)p,(x) (20) 

X 

where x = {xi, . . . ,xn) denotes a defect configuration with defects at sites Xj. The sum 
is restricted to such configurations for which the longest bottleneck has length /* (and 
therefore N > I*). P^{^) is the probability to find the configuration x. 

% There can be more than just one longest bottleneck. 
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Denoting the transport capacity in SB A by Jsba we have 

(J*)(0,L,r) = J^ba(^*) +E'^^'*W^<^W (21) 

X 

with AJ*. (x) = Jf. (x) — Jsba{^*)- The expection value for the corrections to SBA is 
then 

(AJ*)(</),L,r) = ^ Aj;(x)P^(x) = 5]5^'Aj;(x^.)^0(xivO, (22) 

X N* xjv* 

where N* = N — I* is the number of defects besides /* and x^v* denotes the positions of 
these defects. In case of a degeneracy one of the longest bottlenecks is chosen arbitrarily. 

Since P0(x7v) = 0^(1 — 0)'^"^ = (9(0^), the leading correction in O(0) comes from 
configurations with one additional defect: 



(AJ*)(0,L,r) ^ J] Aj;(xi)P^(xi) = (j^Aj;{x,)] P^x,] 



(23) 



where we have used that P^{xi) does not explicitly depend on xi (all allowed defect 
positions are equally probable). 

As long as the longest bottleneck is far from the boundaries, which we can be 
assumed for large systems, the transport capacity does not depend explicitly on its 
position [21] . Hence, instead of Xi we can also use the relative position d of the additional 
defect to the longest bottleneck to characterize the configuration. If the defect is right 
of the longest bottleneck, we have d > 0, else d < 0. Then we obtain the following 
necessary condition for the SBA to work for large systems (L — ^ oo): 

oo ' 

5^Aj;(rf)<oo (24) 

d=—oo 

This condition is fuUfilled if the "bottleneck-bottleneck interaction" AJi*{d) decays 
faster than \d\~^ for large \d\, which is an restriction on the interaction strength of 
defects. In [34J it was shown numerically that in the TASEP this function indeed 
decays faster than |c?|-2, so that (El) is fuUfilled for the TASEP. 

We can further quantify the contribution of the first defect near the longest 
bottleneck as P^{xi) = 0(1 - 0)^"^ = + 0(0^). Since we have to take into account 
defects to the right and the left of the longest bottleneck, we obtain in leading order 



(j*)(0,r 



J^AJ^id) 



(25) 



where contributions with a defect on an adjacent site of the bottleneck (i.e. d = 1 and 
d = —1) do not appear in the sum, since they belong to longer bottlenecks. Note that 
this approximation does not depend on L. 

Unfortunately currently no analytical results for AJf, (rf) are available. Therefore 
we have to rely on the results of MC simulations to test the considerations made in this 
section. We simulated systems with one bottleneck at a position far from the boundaries 
(> 200 sites) and one single defect for several bottleneck lengths /* and defect position d 
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relative to the bottleneck to obtain Jl*{d). The interaction function is then obtained as 
AJ;t ((i) = Ji*{d) — J*,, where Jf* is the transport capacity of a single bottleneck. Since 
AJ^,{d) should decay fast with increasing \d\ (see also [31]), it is sufficient to take into 
account only defects within a finite distance to the bottleneclJ3- In order to obtain the 
expectation value (J*)(0) for arbitrary configurations, one has to average over /* in the 
same manner as in eq. f|T8l) . 

In Fig. [1] we have plotted average values of the transport capacity obtained by MC 
simulations in dependence on the defect density as well for the disordered TASEP and 
the NOSC model. Each data point has been obtained by simulating 50 samples. For 
comparison the results in SBA and the leading order corrections obtained by (|25l) are 
included. We see that while already the SBA appears to be a good approximation, 
the accuracy of the corrections over a wide range of defect densities is astonishing. It 
comes as a surprise that in the TASEP for larger defect densities the leading order 
correction, which takes into account only one additional defect, is extremely accurate. 
This is not expected since for larger there is a higher probability of having more than 
one defect in the vicinity of the longest bottleneck. However, these results indicate that 
the position of other defects beyond the first one do not significantly contribute to the 
transport capacity. Furthermore we see that the deviation of the SBA approaches a 
rather constant value for larger 0, despite the factor in (!25l) . This indicates that for 
larger bottlenecks, the infiuence of single defects on the transport capacity is weaker 
than for small bottlenecks, which is consistent with results in |3lj. For small defect 
densities configurations where the longest bottleneck is near a boundary bottlenecks 
become relevant as was argued in the beginning of the section, thus a deviation of 
the SBA arises in this region, although the distance of other defects from the longest 
bottleneck is large on average. 

6. The phase diagram of disordered driven lattice gases 

The phase diagrams of driven diffusive lattice gases that have exactly one maximum 
all have the same topology. This is based on a maximum current principle and shock 
dynamics |2l Ull |12]. The class of DLGs meeting this condition includes many weakly 
interacting lattice gases, e.g. the TASEP and the NOSC model [361 EZ]. If disorder is 
included, some conceptual problems with the expression phase diagram arise. Usually a 
phase transition is identified by a non-analytic behaviour of a macroscopic quantity. In 
driven lattice gases these can be discontinuities in the density (first order transitions) 
or kinks in the dependence of the current on the system parameters (second order 
transition). Strictly speaking these transitions only occur in infinite systems, since 
non-analyticities can only in the thermodynamic limit. In disordered systems, however, 
there is no unique way of taking the limit L — oo since this can not be done with 
a fixed defect sample and, as we have seen, macroscopic quantities like the transport 
capacity may be sample-dependent. Indeed, the process of taking the thermodynamic 

In our computations we simulated systems from d = —20 to d — 20. 
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Figure 1. First order corrections to SBA as function of the defect density (j> for 
the disordered TASEP (top) and the disordered NOSC model (bottom) using the 
probabihty distribution (|10p for averaging over I* . The slow hopping rates are q = 0.6 
in the TASEP and ujf°^At = 0.32 in the NOSC model. The system size is L = 1000 
in each case. 



limit has to be specified, since it is ambiguous how the "new defect sites" by increasing 
L are included. Enaud et al. [35 , e.g. discussed two possibilities of defining a hmit 
L ^ oo and showed that if this limit is taken by including sites at the boundaries there 
actually is no unique phase transition point if exit rate f3 is fixed and a is varied. For 
infinite systems, according to equation (fTTIl . the length of the maximum bottleneck is 
infinite and thus the transport capacity would be the same as the one of a pure system 
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with hopping rate q, J* = g/4. In this work, however, we are exphcitely interested in 
"finite but large systems" and we are considering ensembles, not individual samples. 
Since the longest bottleneck increases as O(lnL), the transport capacity approaches its 
asymptotic value only logarithmically: J*{L) = q/^ + 0(1/ In L) (see also |23]). For 
finite but large systems we have to take into account terms of the order 0(1/ In L). 
Hence in this view, we want to consider an explicit dependence on the system size and 
cannot take the thermodynamic limit to obtain phase transitions. In [31] it is shown 
that if a single bottleneck is near a boundary, phase separation cannot occur. In this 
case the character of phase transitions is different, since the current is not limited by 
the bottleneck anymore but by the bulk exclusion like in the pure system. In this 
case the phase transition is of second order. On the other hand, if the bottleneck is 
far from the boundaries at a distance d = 0{L) there is not only a sharp kink, but 
also macroscpic phase separation occurs accompanied by a steep increase of the average 
density, indicating a first order transition. 

In Sec. Owe have seen that the average distance of the longest bottleneck from the 
boundaries is 0{L). Hence on average we have a sharp transition for large L. Therefore 
we call this a phase transition for finite but large systems at the critical point ac where 
the current reaches J*, although this point depends on the system size. 

6.1. Effective boundary rates 

Investigations in a TASEP with one and two bottlenecks far from the boundaries 
(distance 0{L)) |Ml[32] showed that the transport capacity only depends on the longer 
bottleneck, while outside of the maximum current phase the current only depends on 
the position of a bottleneck that is near a boundary. This (negative) edge effect is 
considerable for defects not more than ~ 20 sites away from the boundaries. The 
observation of this effect motivates the concept of effective boundary rates: If the system 
is not in the maximum current phase, it can be treated as a pure TASEP with effective 
boundary rates aes,PeH that depend on the distance and size of a bottleneck from the 
boundary and differ from the real boundary rates a and /?. The concept was tested 
in [34] for single bottlenecks and yielded good results. The transition from low- to 
high-density phase was shown to be at the line aes = PeS which in general does not 
correspond to the diagonal a = (3 in the phase diagram. The observations of Enaud et 
al. [35] in the disordered TASEP for different defect samples indicate that the concept 
of effective boundary rates can also be apphed for the disordered TASEP. 

Taking into account defects near the boundary, we can write the current in the 
form J(a) = a{l — a) + AJq,(x). Here a is the entry rate in the low density phase. 
However due to particle-hole symmetrj0 we can transfer this result to (3 and the high 
density phase. The defect configuration {xi,X2, ■■■) is defined in the same manner as in 
Sec. [5l Indeed, taking the expectation value we can proceed analog as in the last section 

* Note that for individual defect samples, particle-hole symmetry is broken, but for large ensembles it 
is restored. 
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to obtain the average corrections in leading order 



(26) 



where di is the position of the first defect and AJa{di) = Ja{di) — a{l — a). Thus the 
corrections by defects near the boundaries are of the same magnitude as the corrections 
to SBA, while the "defect-boundary interaction" AJa{di) is in general not the same as 
the "defect-defect interaction" AJi*{di). Fig. [2] shows that results obtained from 
yield an accurate approximation for the expectation value of the current for low entry 
rates. 

The expectation value of the effective entry rate can then be obtained if the current 
density relation of the pure system J(p) is known. 




Figure 2. Disorder average of the current in dependence of the defect density in the 
disordered TASEP with q — 0.6, a — 0.2 (i.e. J < J*). MC simulations are compared 
with leading order corrections to the pure current obtained by (j26p . 

If the relations aefr(a), PesiP) and their inverses a~^{aes) and P~^{Pes) are known 
as well as the transport capacity J*, we are in principle able to map the problem of 
determining the phase diagram of a disordered system on a pure system with a known 
dependence of the current on the boundary rates J(a,/3): 

(i) If the system current J{aes,Pes) < J* the system globally has the same properties 
as the pure one if one replaces the real boundary rates by the effective ones. 

(ii) At the points in the a— /3-space where J{aes, Pes) = J*, a. phase transition occurs to 
a phase separated phase occurs in which the current is independent on the boundary 
rates and maximal. 
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In particular in the TASEP we can determine the expectation value of the effective 
boundary rates 

(«eff) = \j\-{Jo) (27) 

There is a phase transition from low density to high density phase for aes{oi') = 
PesiP') -v^ a' = «rff (/?eff(/5)) which in general is not on the diagonal a = p. Nonetheless 
we have on average due to particle-hole symmetry (Jq,) = (J^) that leads to a' = P' on 
average. The transition to the phase separated phase is determined by ac{l — etc) = J* 
or Pc{l — Pc) = J*- Unfortunately, we are not able to determine the functions aefr(«), /^efr 
explicitely, since for each a, P we need to obtain a set of functions A Jq,^^ which requires 
much computational effort, as long as no analytical results are available. Nonetheless, 
the concept of effective boundary rates can be used to extract some qualitative properties 
of the phase diagram, though obtaining quantitative results is difficult. 

However, since corrections of the SBA are of the same order as corrections to the 
boundary rates, we can approximate aes ~ a and Pes ~ /9 in order to find ac and Pc- In 
Fig. [3] we plotted the current and the average density in dependence on the entry rate 
a. One observes a steep increase in the average density at the point where the plateau 
begins. Thus we can can characterize this transition as a first order phase transition. 




Figure 3. Mean density (p) and current J in dependence of a for fixed (3 = 0.9, 
q = 0.6 and <j> — 0.05 obtained by MC simulation of a system with L — 2000 and fixed 
defect sample. One observes a steep increase in at the same point where the current 
reaches the plateau. 

Fig. m displays a sketch of the phase diagram of an individual defect sample in the 
disordered TASEP. The transition line between HD and LD is distorted compared to 
the homogeneous case. Taking the disorder average, the transitions are again on the 
diagonal line a = p. The maximum current phase is enlarged and can be characterized 



Driven Lattice Gases with Disorder and Open Boundary Conditions 



18 



as a first order transition compared with the pure system, since a jump of the average 
density occurs at this hne. 
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Figure 4. Schematic phase diagram of the disordered TASEP for a single defect 
sample. It is obtained by taking the phase diagram of the TASEP in dependence of 
tteff and /3eff with one bottleneck that corresponds to the longest one and rescale the 
axis by aeff ~* ol-, /^eff ~* /3- For comparison we included the phase transitions of the 
homogeneous TASEP with p = 1 (dashed lines). 



7. Discussion and Outlook 



We have investigated disorder effects in driven lattice gases using the TASEP with binary 
hopping rates as paradigmatic example. Our results indicate that good approximations 
for the expectation value of the transport capacity (and other quantities) of an ensemble 
of driven diffusive lattice gas systems can be obtained with low effort (no averaging over 
disorder distribution). 

The basic idea stems from the observation that the longest bottleneck (consecutive 
string of slow sites) is the current limiting factor. This suggests the possibility of 
calculating the transport capacity of given defect samples by the Single Bottleneck 
Approximation (SB A). It allows to use known results for systems with only one 
bottleneck, which are usually much better understood than the disordered ones, as an 
efficient and accurate description. With the help of extreme value statistics one obtains 
the probability distribution for the longest bottleneck from which the expectation value 
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of the transport capacity in the SBA can be determined. Since for finite systems only a 
small range of bottleneck lengths gives relevant contributions to the expectation value, 
it is sufficient to have the results for a small number of single bottleneck systems. So 
even for systems for which no analytical results are available, one can get SBA results 
by once simulating a small number of single bottleneck systems. Using these data 
approximations of the transport capacity for arbitrary system size and defect density 
(but fixed transition rates p, q) can be obtained. 

We emphasize that the results obtained here are useful in two different situations: 
a) for a fixed realization of disorder, if the longest bottleneck can be identified, and b) 
for ensembles of systems with fixed density. In the first case, one can directly identify 
the disordered system with the appropriate single-bottleneck case. In the second case, 
which is also relevant for many realistic scenarios, one can use the statistical description 
developed here to obtain predictions for the ensemble. 

The accuracy of the SBA can be systematically improved by taking into account 
various corrections. We found that for small defect densities the most important 
correction is due to the first defect next to the longest bottleneck. It can be expressed in 
terms of functions AJ*(xi) that measure the contribution of a single defect at position 
Xi relative to the bottleneck. A rather general argument indicates that the SBA is 
applicable to a generic driven lattice gases if these functions decay faster than with 
increasing distance. Indeed for both cases studied here explicitly, the disordered TASEP 
and the NOSC model with vanishing Langmuir kinetics, the SBA yields good results 
and thus we can expect it to work even for generic driven lattice gases. Surprisingly, 
the leading order corrections appear to be quite accurate in both systems also for larger 
defect densities (f) ~ 0.5. This indicates that other defects than the first one only 
have very small infiuence on the transport capacity if large bottlenecks are presenlp. 
However for finite systems at small defect densities deviations from the SBA occur that 
cannot be explained by defects near the longest bottleneck. Here interactions of the 
boundaries with (one of) the longest bottleneck can not be neglected which lead to 
relevant deviations (see Sec. [5]). 

We observed that deviations from the current of the pure system also occur if the 
current is less than the transport capacity. However, these deviations are smaller in 
magnitude. This effect is due to defects near the boundaries, which was already shown 
before in systems with single bottlenecks [311 ES]. They can be treated in the same 
manner as corrections to the SBA and we see that in this case the first defect near the 
boundary is the most relevant contribution as expected from the results before. The 
effect can be encompassed in terms of effective boundary rates. In principle for known 
relations between boundary rates and transport capacity, the problem of determining 
the phase diagram can be mapped on a pure system using these quantities instead of 



the ones of the pure systen Jtfl . Though usually it is difficult to determine effective 



The average length of the longest bottleneck increases with increasing defect density, 
ft One still has to be careful since the characteristics of the phases can be different in the disordered 
system, although the topology is the same. 
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boundary rates explicitely that concept is useful to obtain qualitative properties of the 
phase diagram. 

From a theoretical point of view the SBA and its corrections as well as the effective 
boundary rates are interesting since by these concepts disordered systems can be treated 
in terms of systems with single bottlenecks and two-bottleneck systems. These are 
much easier to investigate since one has to consider fixed defect configurations. This 
follows the tradition of statistical physics since microscopic properties of particles as 
well as particle-particle interactions (here "bottleneck-bottleneck interactions" in form 
of the functions AJ(x)) are used to obtain macroscopic quantities using statistics. The 
concept presented in this work is rather generic provided that microscopic properties 
can be obtained. This can for example be done by numerical simulations. 

Driven diffusive systems are used as models for active intracellular transport 
processes. These are characterized by the directed motion of motor proteins on 
microtubules. However, usually the microtubuli are not homogeneous, but there are 
other microtubule associated proteins (MAPs) that are attached to the microtubules 
and can form obstacles that correspond to defects on the modelling level, impeding 
forward movement of motor proteins. One example is the aggregation of tau proteins 
in neurons of organisms suffering of Alzheimer's disease [13|. Furthermore there are 
experiments that show that modified kinesin molecules can immobilize moving kinesins 
which thus form obstacles on the microtubule track For living organisms the current 
of transported objects is a measure for the performance of the transport system, which 
may not fall below a threshold for maintaining cell metabolism and enable cell division. 
Hence, the maximum current is a measure for the transport capacity of a microtubule. 
Since binding and unbinding of molecules to microtubules and kinesin occurs stochasticly 
depending on temperature and concentration, this system meets the criterion of a 
randomly disordered system. In a living organism there can be trillions of microtubules, 
thus the expectation value of the maximum current is a crucial quantity. The defect 
density rather than the individual sample of defects on a microtubule is a measurable 
quantity determined by the concentration of defect molecules and temperature. 

Nonetheless systems with particle conservation in the bulk are not sufficient to 
serve as models for intracellular transport. One crucial property intracellular transport 
exhibits is the attachment and detachment of motor proteins. That means that one has 
to include these effects in the models allowing particle creation and annihilation that 
leads to a spatially varying current. The so called PFF model |5] includes Langmuir 
kinetics to the TASEP and virtually takes into account the attachment and detachment 
of moter proteins. In [38] this model was investigated with one defect site. The NOSC 
model in its original form [36] also includes creation and annihilation of particles and was 
used to model the dynamics of the KIFIA motor protein using some kind of Brownian 
ratchet to perform directed movement. The success of the SBA provoces the assumption 
that the defects locally impose a maximum transport capacity, so that the spatial 
varying current may not exceed the transport capacity at any point. This problem is 
currently under investigation and the results may help to improve our understanding of 
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intracellular transport with particle creation and annihilation in the presence of defects 
or disorder. 

Appendix A. The NOSC model 

The NOSC model is used for modelling the dynamics of KIFIA motor proteins on 
microtubules. These motor proteins can be in a strongly bound state (1) where 
movement parallel to the microtubule is not possible, and a weakly bound state (2) 
where it can diffuse along the microtubule, cyclic transitions between these to states 
leads to a directed net motion using a Brownian ratchet due to an asymmetric binding 
potential. The transitions rules are in the bulk 

Transition to weakly bound state (hydrolysis): 1 — 2 with probability ujh^t 
Transition to strongly bound state on site i: 2 — 1 with probability uj^At 
Transition to strongly bound state on site i + 1: 20 ^ 01 with probability ujfAt (A.l) 
Diffusion to the right: 20 02 with probability UfyAt 

Diffusion to the left: 02 — > 20 with probability UbAt 

at the right boundary (site 1): 

Entry at the left boundary: 0—^1 with probability aAt 
Detachment: 1 — > never (A. 2) 

Diffusion out of the system: 2 — ^> never 

at the left boundary (site L): 

Exit 1^0 with probability (3 At 

Attachment — » 1 never (A. 3) 

Diffusion out of system 2^0 never 

In the original NOSC model [36] also creation and annihilation of particles, i.e. Langmuir 
kinetics, is included. 

Creation: 0-^1 with probability UaAt ,^ ^, 

Annihilation: 1^0 with probability cUdAt ' 

but here we focus on the case Ua = ujd = 0. In our simulations we have considered 
disorder in the forward-rebinding rate uj which is one of the parameters that control 
the average velocity of a particle. The standard parameter values are 0;^^* = 0.145 ms~^ 
for the fast rate and cu^-'"™ = 0.08 ms~^ for the slow rate. In this work we used a timestep 
of At = 4 ms, thus the transition probabilities are obtained by four times the rates. 
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